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ABSTRACT 

Accurate measurements of the mass distribution in galaxy and cluster halos are essential to 
test the cold dark matter (CDM) paradigm. The cosmological model predicts a universal shape 
for the density profile in all halos, independent of halo mass. Its profile has a 'cuspy' centre, 
with no evidence for the constant density core. In this paper we carry out a careful analy- 
sis of twelve galaxy clusters, using Chandra data to compute the mass distribution in each 
system under the assumption of hydrostatic equilibrium. Due to their low concentration, clus- 
ters provide ideal objects for studying the central cusps in dark matter halos. The majority 
of the systems are consistent with the CDM model, but 4 objects exhibit flat inner density 
profiles. We suggest that the flat inner profile found for these clusters is due to an underesti- 
mation of the mass in the cluster centre (rather than any problem with the CDM model), since 
these objects also have a centrally peaked gas mass fraction. We discuss possible causes for 
erroneously low mass measurements in the cores of some systems. 
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1 INTRODUCTION 

In hierarchical collapse models, TV-body simulations predict a uni- 
versal shape for the mass distribution in dark matter halos from 
dwarf galaxies (10^ M0) to massive clusters (10^^ Mq), indepen- 
dent of the value of the cosmological parameters. The density pro- 
file differs strongly from a simple single power-law, predicted in 
early t heoretical work jciunri & G otl 1972;^ Fillmore & Goldreich 
11984 iHoffman & ShahanJ Il985t Iwhite & Zaritsk^^ fl99^ . and 
is centrally concentrated with an inner cusp p(r) oc r~". 
The universality of the structure of dark matter halos formed 
through hierarchical clustering was 'discovered' in simulations per- 
formed bv Navarro et al. (1995, 1996, 1997), hereafter NFW\ 
In these early studies, an inner slope with a — 1.0 pro- 
vided the best-fit to the spherically averaged density p rofiles . 
This was supported by various authors ^ole & Lacey Il996t 
iTormen et alj 1 19971 : iKravtsov et alj Il997 *). Later simulations at 
high er resolution indicated a s t eeper centra l cusp, with 



1.5 jpukushige & Makinj|l997l. I2OOIL l2003l : iMoore et alJll999l: 
iGhigna e t alj200d:lKlvDin et all200lh . A density profile similar to 
the NEW profile, but with a = 1.5 was suggested by Moore et al. 
il999), hereafter M99. 

A thorough understanding of the predictions made by simu- 
lations is essential in order to test the cold dark matter (CDM) 
paradigm against observations. There are several issues which need 
to be addressed by simulators, in particular in the high density cen- 
tral regions of halos where large numbers of particles and fine time 



^ Sim ulations by NFW built on ea rlier, pioneering worlcs by lFrenket alj 
IT988 i. Dubinski & Caribers^ Il99ll) and lCrone etal] ll994l) . which identi- 
fied the absence of a constant density core in dark matter halos. 



resolution are required. In recent work there has been some de- 
bate over the for m of the pro fi le at t h e innermost re s olved radii 
(r < O.Olryj.). [ Power etalJ <2003h . iHavashi et all <2003h and 
iNavarro et alj i20()4ll suggest that the inner density profile does not 
converge to a well-defined power-law, but continues to flatten in- 
wards. Die mand et al- ( 2Q.Q5) . on the other hand, dispute this result, 
continuing to support the existence of a central asymptotic power- 
law to the density profiles. 

A further issue which rem ains unresolved num erically is the 
universality of the core profile. IJing & Sutd <200Cl) . for example, 
find that the inner slope steepens with decreasing halo mass. In 
add ition, using analytica l arguments, Hoffman & Shaham ( 198^ 
and lSver&WhiteHl998h show that the halo profile should depend 
on the power spectrum of initial density fluctuations, with an inner 
slope given by 3(n+3) / (n+4) and 3(n-|-3) / {n+5), respectively, 
where n is the power-law index of the spectrum. 

Nevertheless, there are two robust predictions made by simu- 
lations which may be tested against observations: the density pro- 
file in CDM halos differs strongly from a single power-law, and 
the inner region is cuspy with a power-law slope in the range 
1-0 J; Q ^ 1-5 down to at least 1 per cent of the virial radius. 
In this paper we obtain mass profiles under the assumptions of hy- 
drostatic equilibrium and spherical symmetry for a sample of X-ray 
peaked galaxy clusters using Chandra data. Straightforward depro- 
jection techniques are robust and work well on such clusters. In 
addition, clusters halos have smaller concentrations than galaxies, 
making them ideal objects to study cusps. Work on the inner den- 
sity slope of clusters o f galaxies has also b een carried out in recen t 
stu dies bvlKatavama & Havashida ( 2 004l) . lArabadiis et alJ <2004 
and lPointecouteau et al^ ( t2005») . 



© 0000 RAS 



2 Voigt & Fabian 



A ACDM cosmology with Qa =0.7 and h = 0.7 is adopted. fits|^, although linked between shells. The nominal (measured by 

Ipickev & Lockman 1990) and best-fitting A'^h values are shown in 
Table|2| together with the minimum chi-square of the fit. 



2 CLUSTER SAMPLE 

The cluster sample used is listed in TableQ Observations of objects 
with known short central cooling times, indicating a relatively re- 
laxed object, were obtained from the Chandra archive^. The sam- 
ple of 12 objects was chosen to cover a range in cluster redshift 
(z = 0.02 - 0.46) and ICM temperature (/cbT = 2 - 15keV). 
High redshift clusters observed by Chandra tend to be hotter, and 
in this sense the sample is not statistically complete; however, this 
will not in general affect the analyses carried out, and any biases 
introduced by the limited sample are discussed where relevant. 

Data processing was carried out using the CIAO (Chandra 
Interactive Analysis of Observations) software package available 
from the CXC (Chandra X-ray Centre). Level 2 events files were 
updated to include the latest calibration. Only those X-ray events 
with ASCA grade classifications 0, 2, 3, 4 and 6 were included in 
the cleaned data set. Further screening was carried out to remove 
time intervals contaminated by background flares. The lightcurves 
were analysed using the script LC.CLEAN written by M. Marke- 
vitch and provided by the CXC. The observation length remaining 
after screening (GTI; good time interval) is listed for each object in 
TableQ 



3.3 Deprojection 

A series of tests have been carried out by Ijohnstone et alj i2004h 
to check that pro jet produces the correct results when applied 
to known synthetic data. With the exception of the outermost shell, 
the temperature and density profiles obtained with pro jet match 
very well with the synthetic input profiles. 

An apparent increase in density in the outermost shell re- 
sults from the assumption of zero emission exterior to that region. 
Clearly, there will be some counts attributed to this shell which 
were emitted at larger radii and the emission integral will be over- 
estimated. This effect will be seen out to large radii — even if 
the cluster counts exterior to the outer region are negligible, the 
background emission will contribute to the emission integral — al- 
though will be become less significant as the virial radius is ap- 
proached. 

We note that the measured gas temperature in the penultimate 
annulus is slightly higher when the outer shell is removed from 
the fit, although consistent within the one sigma limits. The tem- 
perature in the outer annulus will be slightly overestimated if the 
temperature rises exterior to this region and underestimated if it 
falls. 



3 SPECTRAL ANALYSIS 
3.1 Extracting spectra 

For the spectral analysis, each cluster was divided into circular an- 
nuli around the X-ray emission peak, with strong point sources 
identified by eye and masked out. Regions disturbed by, for ex- 
ample, merging subclusters and radio lobes, were also removed 
(see Section 3.4). A spectrum was extracted for each annulus us- 
ing the CIAO DMEXTRACT tool. The spectra were binned to con- 
tain at least 20 counts per PHA channel to enable the use of chi- 
square statistics. Ancillary-response and response matrices were 
constructed using the CIAO MKWARF and MKRMF programs. 



3.4 Masked regions 

To obtain an accurate measure of the cluster mass distribution we 
are interested in the ambient gas properties. For the spectral anal- 
yses, the following features were therefore masked out from the 
images: a shock front in 2A 0335+096 spanning a sector between 
150 — 210° 50 — 60 arcsec from the cluster centre; radio lobes in 
Abell 478; the southeast quadrant in RXJ 1347.5—1145 containing 
a bright subclump; the bright filament extending southwards from 
the core in Abell 1795. 



3.2 The spectral model 

Spectra were fi t in the 0.5 — 7.0 keV energy range using the XSPEC 
jAmau?199^ software package. The emission from each shell was 
modelled using the mekal iMeweetal. 198 ^ plasma emissio n 
code, incorporating the Fe L c alculations of iLi edahl et alj (l99^, 
and a bsorbed by the phabs iBalucinska-Church & McCammonI 
Il992h photoelectric absorption code to take into account galactic 
absorption along the line-of-sight. The emission was deprojected 
using the pro j ct model provided in XSPEC. 

The free parameters in the model were the temperature, metal- 
licity and emission measure. The elements were assumed to be 
present in the solar ratios measured by Anders & G revesse) il989h 
and the abundance allowed to vary between shells. The Galac- 
tic absorption column density was left as a free parameter in the 



^ Observations using the ACIS-S3 detector only were used. 



3.5 Spectral fits 

The reduced chi-square of the fits obtained using the model 
pro jct*phabs(mekal) are less than ~1.1 for the majority of 
objects in the sample (see Table|2j, showing that the X-ray emission 
is well described by a single phase plasma at each radius. We note 
that the reduced chi-square is greatest for the fits to 2A 0335+096 
and Abell 478. 



The low energy quantum efficiency of the ACIS chips has been contin- 
uously degraded since launch. Tools have been provided by the CXC to 
account for the loss in effective area, although whether or not the correc- 
tion is applied has a negligible effect on the temperature and emission in- 
tegral profiles when the galactic absorption column density is left as a free 
parameter in the fits. With the correction applied, the best-fitting A^h is 
often much less than the nominal value, suggesting t hat the tools ten d to 
'over-correct' the data. Following the work of Voigt & FabianI l2004l) and 
j^irzan et al. 1 2()o4, we use uncorrected data in order to avoid best-fitting 
galactic absorption measurements consistent with zero. 
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Cluster 


Red.shift 


Obs. date 


Exps. time 


GTI 


Emission peak (J2000) 




Da 








(ks) 


(ks) 


I^A 


Dec 


(Mpc) 


(kpc arcsec"^ ) 


Abell 3112 


0.0750 


2001 Sep 15 


17.1 


13.5 


03 17 57.7 


-44 14 16.9 


335.2 


1.4 


2A 0335+096 


0.0347 


2000 Sep 06 


20.0 


18.1 


03 38 40.6 


+09 58 11.0 


140.6 


0.7 


Abell 478 


0.0880 


2001 Jan 27 


42.9 


38.9 


04 13 25.2 


+ 10 27 53.9 


396.9 


1.6 


PKS 0745-191 


0.1028 


2001 Jun 16 


17.9 


14.6 


07 47 31.2 


-19 17 38.8 


468.5 


1.9 


RXJ 1347.5-1145(1) 


0.4510 


2000 Mar 5 


9.1 


8.0 


13 47 30.5 


-11 45 09.5 


2492.9 


5.7 


RXJ 1347.5-1145(2) 


0.4510 


2000 Apr 29 


10.1 


7.2 


13 47 30.5 


-11 45 09.5 


2492.9 


5.7 


Abell 1795 


0.0632 


2000 Mar 21 


19.7 


15.6 


Lj 4o JZ..J 


1 1^ 11 Q 


280.0 


1.2 


ADell IfSjJ 


U.ZjzJ 


1 QOQ 1 1 

1 yyy uec 1 1 






14 01 01.9 


+02 52 43.4 


1 TAT A 


1 a 

3.y 


Abell 3581 


0.0218 


2001 Jun 07 


7.3 


6.0 


14 07 29.8 


-27 01 04.2 


93.6 


0.4 


Abell 2029 


0.0767 


2000 Apr 12 


19.9 


19.8 


15 10 56.1 


+05 44 40.6 


343.2 


1.1 


RXJ 1532.9+3021 


0.3615 


2001 Aug 26 


9.5 


6.2 


15 32 53.8 


+30 20 58.5 


1916.3 


5.0 


MS 2137.3-2353 


0.3130 


2000 Dec 10 


44.2 


19.5 


21 40 15.2 


-23 39 39.9 


1618.8 


4.7 


Sersic 159-03 (AS 1101) 


0.0564 


2001 Aug 13 


10.1 


9.8 


23 13 58.3 


-42 43 35.0 


248.7 


1.1 



Table 1. Summary of the Chandra observations. Cluster redshift, observation date, exposure time, good time interval. X-ray emission peak, luminosity distance 
and angular scale. Both observations available for RXJ 1347.5— 1 145 were used. 



Cluster 


Fitted A^H 
(lO^O cm-2) 


Galactic A^h 
(10^0 cm-2) 




xl 


Abell 3112 


5 42+"'22 
■^■^■^-0.22 


1.95 


1482(1388) 


1.07 


2A 0335+096 


11 +0- 18 
^'■^^-0.18 


17.8 


2682(1968) 


1.36 


Abell 478 


33 4^+0.16 

JJ'+D-o.ie 


15.2 


4243(3458) 


1.23 


PKS 0745-191 


47 77+0-35 


42.4 


2486(2251) 


1.10 


RXJ 1347.5-1145 


;r 1O + 0.72 
31S_0.41 


4.85 


805(859) 


0.94 


Abell 1795 


, oT+0.13 
^■^^-0.13 


1.19 


2110(1861) 


1.13 


Abell 1835 


2./»_o.24 


2.32 


1551(1445) 


1.07 


Abell 3581 


n g7+0.53 
^■°'-0.50 


4.52 


716(653) 


1.10 


Abell 2029 




3.05 


2878(2555) 


1.13 


RXJ 1532.9+3021 


6 14+0-68 


2.16 


403(379) 


1.06 


MS 2137.3-2353 


4 63+0.44 


3.55 


596(541) 


1.10 


Sersic 159-03 


gQ7+0.39 
°'^"-0.39 


1.79 


917(844) 


1.09 



Table 2. Best-fitting absorption column density, A'^pj, chi- square and reduced chi-squ are for the model pro jct*phab3(mekal). The nominal galactic 
absoiption column along the line-of-sight to the cluster from lDickev & LockmanI ll990t) is also tabulated. 



4 TEMPERATURE, DENSITY AND PRESSURE 
PROFILES 

Gas temperature, density and pressure profiles are shown for a rep- 
resentative sample of clusters in FigQ 

4.1 Temperature and density 

The temperature distribution is obtained directly from the spectral 
fits and the density profile from the normalization of the mekal 
spectrum, K = _E//(4 x 10^^ttDI{1 + zf), where EI{= 
J neUHdV) is the emission integral. The electron number density 
at the centre of each shell is estimated as (rie) = {1.2EI /Vy^^ , 
where we have used the relation rie ~ 1.2nH and V is the volume 
of the shell. 

4.2 Pressure 

At the high temperatures and low densities in the ICM the ap- 
propriate pressure equation of state is the ideal gas law, given by 
P = nksT, where n ~ 1.92nc is the number of particles per unit 



volume. The electron pressure, Pc, at the midpoint of each shell 
is computed by multiplying the electron density by the tempera- 
ture. The uncertainties in Tie and T are added in quadrature to find 
the uncertainty in Pc (i-e. the density and temperature are assumed 
to be independent. We note that contour plots of shell normaliza- 
tion against temperature are approximately circular, showing the 
assumption is a reasonable one). 

5 MASS MEASUREMENTS 
5.1 Total mass 

Determining masses using X-ray data allows the cluster mass to 
be computed as a function of radius. The gas is assumed to be in 
hydrostatic equilibrium within the cluster potential. The equation 
of hydrostatic equilibrium for a spherically symmetric system is 
given by 

G pgas,i \dr J . 
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Figure 1. Deprojected gas temperature, density and pressure profiles for a representative sample of clusters. Projected gas temperature profiles are shown in 
green. 



where P is the thermal pressure. 



5.1.1 Calculation 



The total mass enclosed within radius r is calculated from Equa- 
tionQ where r-j is the radial distance from the centre of the cluster 
to the midpoint between two consecutive shells. The pressure gra- 
dient at the midpoint is estimated as 



fdP 



(2) 



where Pi is the pressure at the centre of the ith shell, Vi is the radial 
distance to the centre of the ith shell and rj — (rt + ri+i)/2. The 
gas density at rj is calculated using linear interpolation, such that 



Pgas(ri + l) +Pgas(r-«) 



(3) 



where Pgas.i, the gas density at the centre of the ith shell, is 
found from the electron number density using the relation Pgas,i ~ 
^mHWgas,* ~ 1.92/imHnc,i. The uncertainties are calculated us- 
ing a Monte Carlo technique whereby both the pressure and density 
data are perturbed 1000 times. PKS 0745 — 191 does not exhibit a 
smooth pressure profile in the centre and so measurements made 
using Equation|2|do not produce a monotonically increasing mass 
profile. We compute the mass profile in the inner region of this ob- 
ject by binning together the first and second and the third and fourth 
data points in the pressure and gas density profiles. 

Computing the mass profile involves determining the pressure 
gradient at the midpoint between the centres of consecutive shells. 
By estimating the gradient using Equation |2| we effectively fit a 
straight line between adjacent data points. (The accuracy of the 
method may therefore be improved by increasing the number of 
data points). Neighbouring mass points may be slightly correlated, 
but data points further away than one will not be correlated with 
one another. Several authors fit a parametric model to the temper- 
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ature and gas density profiles and use this to compute the mass 
profile. The data points in this case are not independent and finding 
the best-fit model to the mass profile using chi-square statistics will 
be statistically invalid. The method used here, made possible due 
to the high spatial resolution of the Chandra observatory, therefore 
provides an improvement over previous analyses. 

The mass profiles are shown for the sample in Fig. |2| 
The outer mass data point is removed in all but three clusters 
(RXJ 1347.5-1145, RXJ 1532.9+3021 and MS 2137.3-2353). 
This is because the artificially high gas density and pressure in the 
outer bin'' causes the outer mass data point to be significantly un- 
derestimated in the majority of objects. For the high redshift sys- 
tems (2 > 0.3) the outer mass data point is estimated to be accu- 
rate to within a few per cent. The accuracy of the outer data point 
is assessed by calculating the mass profile from the data obtained 
both with and without the outer shell included in the spectral anal- 
ysis. As an example, the mass profiles obtained for Abell 1795 and 
RXJ 1347.5— 1 145 are shown in Fig.|3| The percentage decrease in 
the mass measured in the penultimate shell is less than 1 per cent 
for RXJ 1347.5-1 145 and about 50 per cent for Abell 1795 when 
the outer annulus is not included in the spectral fit. 

Removing the outer mass data point in clusters such as 
Abell 1795 before analysing the mass profile is vital. Current struc- 
ture formation models predict that the mass profile should 'turn- 
over' at some characteristic radius. If an erroneously low mass data 
point is included in the outer region then this will wrongly indicate 
that the profile is curving away from a power-law. 



5.1.2 tensing masses 

Mass measurements from strong lensing analyses are also plotted 
in Fig.|2|for PKS 0745-191, R XJ 1347.4-1145, Abell 1835 and 
MS 2137.3-2353 <AlleiJll998h . In each case the lensing mass is 
larger than the X-ray mass by a fa ctor of between 1.5—3. This 
discrepancy is also found by Allen ( 1998) if a single phase spec- 
tral model (i.e. no cooling flow) is adopted. It is difficult at this 
stage to draw any strong conclusions about the difference in the 
measurements. The lensing results are quoted without uncertain- 
ties and there are several difficulties with this method, such as 
adopting the correct geometry, obtaining accurate arc redshifts and 
the unknown presence of secondary matter along the line of sight 
iWambsganss et all2005h . Similarly, the masses obtained from the 
X-ray analysis are also dependent on several assumptions, includ- 
ing the spectral and geometrical models adopted. It will be impor- 
tant in future work to improve on both methods until an agreement 
is reached. 



5.2 Gas mass and gas mass fraction 

The integrated gas mass profile is calculated using the sum 

^ 4 

Mgas(< rfc) = 3'^(''fc - '^Li)pgaB('-i), (4) 

i = l 

where is the radial distance to the outer boundary of the ith shell. 
^ See Section 3.3. 
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The gas mass fraction is the ratio of the total gas mass to the 
total gravitating mass within a fixed volume 



/gas — 



« r,) 



Mtot(< Tj) 



(5) 



The total gas mass within rj is computed from the gas mass profiles 
using linear interpolation. The gas mass fraction profiles are plotted 
in Fig.|2| 

We note that the total integrated mass within a particular vol- 
ume, calculated from the equation of hydrostatic equilibrium, is 
dependent upon the pressure gradient and gas density at that radius 
only, and is completely unaffected by the regions interior (or exte- 
rior) to that radius. For the gas mass profile, on the other hand, the 
integrated mass is calculated by summing from the centre outwards 
and any error in the measurement at small radii will propagate out 
to larger radii. However, the gas mass at small radii is much less 
than at large radii and any uncertainty in the measurements in the 
core are unlikely to have a significant effect on the gas mass profile 
further out. 



6 DARK MATTER PROFILES FROM SIMULATIONS 

The generalized equation for the density distribution in dark matter 
halos, including all models with a 'cuspy' core, is given by 



p{x) = 



Po 



(1 + xT) 



(/3-c)/7' 



(6) 



where x = r/rg. The scale radius, r-g, is a free parameter. The 
indices (a,/3,7) for the NFW and M99 models are (1.0,3.0,1.0) and 
(1.5,3.0,1.5), respectively. A simple analytic model for the mass is 
easily obtained for each of these special cases, with 

{4,TTpori\\n{l + x) - for NFW 
M(< x)^{ f"-' ^ [ ^ ^ ' i+^J (7) 

Ivrpors ln(l + x^'^) for M99 

A generalized mass profile, valid for any values of the indices, 
may be written 



3 a;"^ 2-^1(03/7, (/3 - a)h; (03 + l)h\ -^'') 



M{< x) = 47rpor, 

L13 

(8) 

where 03 = 3 — a. The hypergeometric function 2F\ may be eval- 
uated numerically. We refer to this model as M2F1. 



6.1 Definition of concentration 

The mass enclosed within a region r a may be written in terms of 
the critical density of the Universe, 



(9) 



where Zf is the cluster formation redshift, assumed equal to the 
observed redshift, such that 



Ma 



4 3 

-Trr-APcritA, 



(10) 



where A is the density contrast. The Hubble parameter, H{z) = 
lQQhE{z) kms~^ Mpc^^ , is derived from the Friedmann equa- 
tion. In a flat Universe E{z)'^ = r2m(l + z)^ + ^a. 



Voigt & Fabian 
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500 1000 




Figure 2. Total (dark plus baryonic) mass profiles (upper) and gas mass fraction profiles (lower). The error bars show the 1 sigma uncertainties in the 
measurements. The solid hne shows the best-fit M2-F1 model to the data, with a a free pai'ameter in the fit. The dashed line shows the best-fit NFW model 
(a = 1) to the data. For the latter model the inner 30 kpc is removed from the fit to the profiles of 2A 0335-f 096, Abell 478, PKS 0745-191 and Sersic 159-03 
(see text). The best-fit scale radius and 1 sigma limits, obtained with a = 1, are shown by solid and dotted vertical lines respectively. The arrow indicates the 
best-fi t scale radius for the M2F1 model. (In the cases where an arrow is not shown the best-fit model tends to a power-law). Lensing masses taken from .Aliei] 
Il998t) are indicated (green filled squares). 
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Figure 3. Mass profiles computed using data from spectral fits obtained with (green filled circles) and without (black open circles) the outer shell included. 



The concentration, ca, is defined such that ca ~ ta/ts. Us- 
ing this definition the normalization of the density profiles, po, may 
be written as the product of the critical density and a dimensionless 
parameter, 5c- Substituting xa = rA/rs ~ ca into the mass for- 
mulae above (Equations|7|and |8j at r = ta and setting them equal 
to EQuation I lOl vields the following expressions for Sc 



3 ln(l+CA)-CA/{l + CA) 
A 

2 ln(l + 4/^) 

A eA"3 



3 2Fi{a3/j,{l3-a)/j;(a3+'y)/f,-cl) 



for NFW 
for M99 

forM2-Fi 



(11) 



6.2 Definition of the virial radius 

The virial radius separates the region where the cluster is in hydro- 
static equilibrium from where matter is still infalli ng. For Qm + 
= 1, rA = rvir ifA = ^vir = 178f7°;''^ jLacev & Cold 
Il993l lEke et ai]|l996l Il99l) . For Qn, = 0.3, Avir = 104. In an 
Einstein-de Sitter Universe, Qm = 1.0 and Avir ~ 178. A density 
contrast of A = 200 (or A = 178) is still often used in the literature 
and Tvir and r2oo (or rns) written interchangeably. For purposes 
of comparison with other work in the literature we adopt A = 200. 



7 MODELLING THE OBSERVED MASS PROFILES 

CDM simulations make two robust predictions concerning the dis- 
tribution of mass in dark matter halos: the shape of the mass profile 
differs strongly from a power-law, and the power-law slope of the 
inner d ensity profile lies in the range 1.0 < a < 1.5 (see, for ex- 
ample, |Nw^^^]|20o3)- These properties are investigated for 
the observed mass profiles below. 



7.1 General shape of the mass profile 

The mass profiles in Fig.|2|are fit with the generalized mass model 
in Equation|8| with (3 and 7 fixed at 1.0. Fig.|4|shows the variation 
in concentration parameter, scale radius, virial radius (r2oo) and 
chi-square with inner slope as a is stepped through values between 
0.0 and 2.0 and the fit is minimized with respect to c and r^. (Note 
that the scale radius and concentration are highly negatively corre- 
lated). The virial radius is computed by substituting ca = ta/tb 
into Eguation ll II 



If the best-fit concentration is less than 1.0 for a particular 
value of a, then the mass model does not agree with CDM pre- 
dictions i.e. the characteristic radius, rs, at which the profile 'turns- 
over' is greater than the virial radius. In Table|3|we show the best-fit 
parameters for the M2F\ and power-law (A/pow oc r^~") models. 
When the best-fit concentration for the M2F1 model is c ^ 1.0 
the model tends to a power-law. We refer to the best-fit model as 
a power-law when the best-fit concentration is less then 1.0. Al- 
though not strictly true, the description serves as a useful distinc- 
tion between profiles which are consistent (c > 1) and those which 
are inconsistent (c < 1) with CDM model predictions. 

With a a free parameter in the fits, the concentration is con- 
sistent with both c > 1 and c < 1 within the one sigma uncertain- 
ties for half the objects analysed (Abell 3112, RXJ 1347.5-1145, 
Abell 1795, Abell 3581, RXJ 1532.9+3021 and MS 2137.3-2353; 
see shaded regions in Fig. |4}. In four clusters (2A 0335+096, 
Abell 478, PKS 0745-191 and Sersic 159-03) the best-fit con- 
centration is greater than 1.0 at the 68 per cent confidence level; 
although the best-fit inner slope in these objects is shallower than 
predicted by CDM simulations. Abell 1835 and Abell 2029 are the 
only two objects for which a power-law is ruled out at the la level 
and the inner slope is consistent with the CDM model. 



7.2 Inner slope measurements 

The best-fit inner slope and la uncertainties are plotted for each 
object in Fig.|5| Also shown are the 3a upper limits on the slope. 
Clusters with a best-fit concentration less than 1.0 are indicated. 
It is clear from the plot that 8 out of the 12 clusters analysed are 
consistent with the CDM model. Furthermore, the points appear to 
scatter around a = 1.0, rather than a — 1.5. An inner logarithmic 
slope of a = 2 is ruled out at the 3a level in all objects except 
Abell 3581; the mass profiles would therefore seem to rule out the 
single isothermal sp here (SIS) model which ha s been widely used 
in the literature (e.g. lBinnev & Tremaine^l987^ . We note that there 
is no obvious correlation between the cluster mass and the best-fit 
inner slope, although the uncertainties are too large to justify any 
correlation tests on the data. 

The four clusters which are inconsistent with an inner 
slope of 1.0 at the la level are 2A 0335+096, Abell 478, 
PKS 0745-191 and Sersic 159-03, with Abell 478, Ser- 
sic 159-03 and PKS 0745-191 consistent with a = 1.0 at the 
2 sigma level. Only 2A 0335+096 disagrees with CDM predic- 
tions with 99.7 per cent confidence. We note that this is also the 
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object for which the spectral model used provides the poorest fit 
(see Section 3.5) and, as discussed below, the temperature distribu- 
tion in the centre is highly asymmetric between the northern and 
southern sectors of the cluster. The presence of an asymmetric tem- 
pe rature distribution in 2 A 0335+096 has been shown previously 
bv lMazzotta et alJj2003h . 

II I \ 
Katavam a & Hayashida ( 2004) find that six 

out of the twenty clusters in their sample (including Abell 478, 
PKS 0745-191 and 2A 0335+096) have a slope lower than unity 
at the 90 per cent confidence level. The authors interpret this result 
as being in agreement with the 'flat core problem' found in obser- 
vations of low mass objects. 

The method used for computing the mass profiles relies on 
the assumptions of hydrostatic equilibrium and spherical symme- 
try. Whether or not the former condition is met is difficult to as- 
certain. Katavama & Havashida (2004) touch on this problem by 
highlighting objects which exhibit obvious structure in their cores, 
suggesting that this may indicate a break from hydrostatic equilib- 



rium, as well as the possibility of affecting the ambient temperature 
and density profile measurements. They conclude that some objects 
have flat inner slopes (Abell 478, PKS 0745-191, ZW 3146), even 
when those with central structure (2A 0335+096 and Abell 2597) 
are removed. 

In Fig|5|we plot the best-fit inner density slopes when the in- 
ner ~ 30 kpc is removed from the fits to 2A 0335+096, Abell 478, 
PKS 0745-191 and Sersic 159-03. It is clear that it is the inner- 
most data point in each of these clusters which is forcing the best-fit 
to be flat. The question is 'how far can we trust data from the cen- 
tral 30 kpc?'. The challenge will be to determine the extent of any 
non-thermal pressure support in the central regions of these objects. 

7.3 Symmetry of the temperature distribution 

One clue that the gas is disturbed and not in hydrostatic equilibrium 
may be suggested by asymmetries in the temperature distribution. 
Here we investigate whether there is any correlation between the 
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Figure 4. Plots showing variation in concentration, scale radius, virial radius and Ax^ (= — Xmin' where Xmin minimum chi-square) for values 
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cent) confidence levels for alpha. The grey shaded regions show the parameter values which lie within the one-sigma confidence interval. Note that when the 
concentration drops below 1.0 the scale radius is larger than the virial radius and the profile tends to a power-law. Fit parameters determined with the inner 30 
kpc excluded are shown with dashed lines, and the one-sigma confidence limits shown by green dashed lines. (If only one green dashed line is shown then the 
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Figure 5. Left: best-fit inner slope. Right: best-fit inner slope with inner 30 kpc removed in 2A 0335+096, Abell 478, PKS 0745-191 and Sersic 159-03. 
The error bars show the 1 sigma uncertainties and the downward arrows the 3 sigma upper hmits. An open circle is used to indicate clusters with a best-fit 
concentration less than 1.0. The clusters are plotted from left to right in the same order as they are listed (from top to bottom) in TablefTI 



degree of symmetry in the temperature distribution and the inner 
slope of the mass profile. 

Fig|6|shows the temperature distributions obtained for each of 
the clusters showing a flat inner density profile, together with two 
clusters which are consistent with the CDM mass model. 

The temperature distribution of 2A 0335+096 is clearly asym- 
metric within the central ~50 kpc. This is unsurprising given that 
[Mazzotta et al. ( 2003) have shown the central region contains many 
blobs of gas at different projected temperatures. The asymmetry 
does not, therefore, necessarily mean that the gas is not in hydro- 



static equilibrium, but does suggest that the temperature distribu- 
tion obtained from the spectral analysis may not be that of the 
ambient gas. In this case, the ce ntral mass measuremen ts obtained 
may be incorrect. We note that lMazzotta et all i2003h find a sig- 
nificant asymmetry in the temperature distribution at radii > 100 
kpc. We also find that the southern half of the cluster is hotter than 
the northern half in the outer region, but the effect is not so pro- 
nounced. There are sev eral differences betwee n the analysis carried 
out here and the one bv lMazzottaetalJ<200 j> which may give rise 
to a discrepancy in the measured temperature distributions. Firstly, 



© 0000 RAS, MNRAS 000, 000-000 



12 Voigt & Fabian 



Cluster 




a 


C200 


\A2F1 model 
rs 
(Mpc) 


''200 
(Mpc) 


X^(dof) 


Power-law model 


Abell 3112 




1 T 1 +0.15 
^■^ -0.55 






4.6 


0.8(3) 


1 Ti+0.15 
1 -'1-0.16 


0.8(4) 


2A 0335+096 


QQQ + 0.15 


12 
15 


0.14 
0.08 


1.7 
1.2 


6.5(2) 
2.3(1) 


55+"°? 

~^ — 0.09 
■ 05+0.21 
l-"-'-0.19 


12.3(3) 
0.7(2) 


Abell 478 




40+0.45 
11^-0.63 


9.0 


0.22 


2.0 
9.2 


6.7(5) 
5.4(4) 


04+0 10 
"■^^-0.10 

11^-0. 16 


7.8(6) 
5 4('5') 


PKS 0745- 


191 


U.UU p pp 
U.UU p pp 


14 
18 


0.13 
0.09 


1.8 
1.6 


2.5(2) 
0.4(1) 


091+0.15 

"■'^ -0.16 

, 09+0.22 

l-'^-0.22 


6 in^ 

0.9(2) 


RXJ 1347.5 


-1145 


1 W+°'^^ 
1-1"-0.41 


3.2 


0.90 


2.9 


3.0(1) 


i-'"-o.io 


3.4(2) 


Abell 1795 




1 24+0'i5 
1-^^-0.56 


— 


— 


5.8 


2.6(2) 


1 24+015 
1-^^-0.16 


2.6(3) 


Abell 1835 




0.52t»;g 


9.2 


0.20 


1.8 


5.9(3) 


1-1^-0.11 


7.5(4) 


Abell 3581 




0.88l»;65 


12 


0.05 


0.7 


1.3(3) 


1 34+0.22 
1-^^-0.24 


1.5(4) 


Abell 2029 




(,74+0.34 
"■'^-0.61 


8.9 


0.21 


1.9 


2.4(4) 


, ,n+0.09 
ll-'-0.09 


A 1 

4.ip) 


RXJ 1532.9 


+3021 


O5; + 0.49 


4.3 


0.46 


2.0 


1.2(2) 


, ,0+0.20 

i-i'^-o.w 


1.6(3) 


MS 2137.3- 


-2352 


1 25+'''l* 
'■'^-'-0.77 






5.9 


0.2(1) 


, 25+0.14 
l-^-'-0.15 


0.2(2) 


Sersic 159- 


03 


QOO+075 


17 
3.4 


0.06 
0.42 


1.1 
1.4 


3.4(3) 
2.5(2) 


,09+0.14 

1"^-0.15 
, ,0+0.21 
l.jb_Q 21 


6.0(4) 
2.5(3) 



Table 3. Best-fit slope a and corresponding chi-square for the M2F1 and power-law models. The second row for 2A 0335+096, Abell 478, PKS 0745— 191 
and Sersic 159—03 shows the best-fit with the inner 30 kpc excluded. The best-fit concentration, scale radius and virial radius are shown for the M2F1 model. 
The concentration and scale radius are not shown for clusters where the best-fit tends to a power-law (i.e. c << 1.0). 



iMazzotta e t 31.72003*) extract spectra in elliptical rather than circu- 
lar armuli and are therefore probing different regions of the cluster. 
If the temperature distribution is not symmetrical then this may lead 
to significant differences in the temperatures measured. In addition, 
[Mazzotta et al. ( 2003) plot the projected temperatures and allow the 
absorption column density to vary between shells. 

For the rest of the sample there is little evidence for a corre- 
lation between the inner slope and the degree of symmetry in the 
temperature distribution. 

7.4 Observed gas mass fraction profiles 

The gas mass fraction profiles are shown in Fig. |2| The pro- 
files were calculated using the relation /gaa(< r) = A/gas (< 
r)/Mtot(< r), where A/tot is the total model-independent mass. 
The profiles are plotted for the whole sample in Fig|7| with the x- 
axis scaled by r2 x lo* (see Section 7.4). The profiles show a variety 
of behaviour within r ^ 0.02rvir, with an apparent flattening at 
/gas ~ 0.1 outside this region, similar to the results of lAllen et alJ 

7.4.1 Inner region 

A high value of /gas, relative to the values outside r ^ 0.02rvir, is 
seen in the central bin of 2A 0335+096, Abell 478, PKS 0745-191 
and Sersic 159—03. These are the only objects in the sample 
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Figure 7. Gas mass fraction profiles for the sample with the x-axis scaled 
by r2 X iQi (see TableBl. The colour-coding is the same as used in Figlsl 

for which the best-fit inner slope, a, to the mass profile is less 
than 1.0 (see Fig|5j, suggesting a link between a high measured 
central gas mass fraction and an apparently shallow inner mass 
density profile. The central rise in /gas is clearly seen in Fig. |8| 
where we plot scaled gas mass fraction profiles for the three 
clusters where the /gas peak is most obvious (2A 0335+096, 
Abell 478 and PKS 0745-191), together with two low redshift 
clusters (Abell 3112 and Abell 1795), whose best-fit mass pro- 
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Figure 6. Plots showing the temperature distribution extracted from the northern (filled circles) and southern (open squares) sectors of the cluster. Also shown 
is the temperature obtained from the whole annulus at each radius (green). 



files are consistent with a = 1.0, for comparison. The latter two 

objects show an increase in fgas, from the centre outwards, as ex- 

I 1 

pected from CDM simulations <Eke et al. 1998). With the excep- 
tion of Abell 3581, whose profile drops towards the centre, the re- 
maining objects in the sample have relatively flat /gas profiles in 
the inner region. 



iKatavama & Havashidal i2004) find a negative correlation be- 
tween a and the value of /gas near the centre of the cluster They 
interpret this as a tendency for clusters with gas-rich cores to have 
flat central density profiles. It seems just as likely that the high gas 



mass fraction measured in the central parts of these objects results 
from an underestimation of the total mass^, supporting the presence 
of non-thermal pressure in the core or suggesting incorrect spectral 
modelling of the data (note that the spectral fit to the emission from 
the centre of 2A 0335+096 is poor). 



^ If the true gas mass fraction profiles of the objects studied here with best- 
fit o < 1.0 are either flat or drop into the centre then this implies that the 
total mass in the inner region has been underestimated by at least a factor 
of 2. The total mass profiles would then be consistent with a ~ 1.0. 
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Figure 8. Gas mass fraction profiles for Abell 3112 (black), 2A 0335+096 
(red), Abell 478 (green), PKS 0745-191 (blue) and Abell 1795 (yellow). 



7.5 NFW model fits 

The majority of tlie clusters in tfie sample are consistent with the 
NFW mass model i.e. a = 1.0 (see Fig.|5}- Here we find the best-fit 
scale radius and concentration for each cluster assuming an NFW 
profile. (For clusters where the NFW model is excluded at the 1 
sigma level, we remove the inner 30 kpc from the fit). It is important 
to stress that attempting to measure is pushing the data to its 
limits since the best-fit scale radius in general lies beyond the radius 
at which mass measurements can be made (see Fig.|2j- The values 
obtained for Vs are therefore strongly dependent on the outermost 
datapoint in the mass profile. It also makes the value difficult to 
constrain and, while some of the clusters show evidence for a turn- 
over in the mass profile, more precise determination of the scale 
radius would require mass measurements out to a larger fraction of 
the virial radius. 

We compute the concentration at three overdensities: A — 
2 X 10*, 2500 and 200. For A = 2 x 10* the value obtained for 
rA is independent of the accuracy of Vs since in each case ta lies 
within the radius at which mass measurements have been obtained. 
However, for A = 200 the value obtained for r2oo is highly de- 
pendent on the accuracy of rs (since it involves extrapolating the 
NFW model using rs). The values for r2oo are shown for inter- 
est, but are not used for any subsequent analysis. The concentration 
is also computed at A = 2500. This allows comparison of our 
results for PKS 0745-191, Abell 183 5, RXJ 1347 5-11 45 and 
MS 2137.3-2352 with those found bv lAllen et ail J200ll ). (Note 
that for the more massive clusters r25oo lies within the radius at 
which mass measurements have been obtained). The data used here 
for PKS 0745-1 91, Abell 1835 an d RXJ 1347.5-1145 are the 
same data used bv lAUen et alji200ll) . The values of r25oo are con- 
sistent between the two studies within the 1 sigma uncertainties. 
However, Allen et al. (2001) generally find a higher concentration 
(and therefore lower scale radius) than the results found here. The 
difference in values is not surprising given the discussion above 
concerning the difficulty of measuring rg. 

The best-fit parameters are shown in Table |4| The mass en- 
closed within a density contrast of 2 x 10* is computed using Equa- 
tion^] The cluster mass profiles, with the a;-axis scaled by ta, 
are shown in Fig.|9| The j/-axis for each mass profile is scaled by 
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Figure 9. Scaled mass profiles: total (filled circles) and gas (open circles) 
mass profiles. The colour-coding is the same as used in Figlsl 



Merit = (4/3)7rrAPcrit A. Each mass profile should therefore co- 
incide if the concentration varies little from halo to halo®. 

A'^-body simulations predict a weak dependence of halo con- 
centration on mass, with Cvir oc M^" and a ~ 0.1 for all ha- 
los at a given redshift. In Fig.^|we plot Cj^io* against M2y^iQ4. 
(Note that we do not plot C200 against M200 here since the calcula- 
tion of A/200 depends on r2oo which is itself an uncertain num- 
ber). The concentration is scaled to a redshift of zero using the 
concentrati on-redshift d ependence found in A'^-body simulations 
by Bulloc k et alJ i200lh . The relation predicts that for haloes of 
the same mass the concentration is proportional to {1 + z)~^ . We 
show the best-fit pow er-law model to the data along with the the- 
oretical prediction of Bullo ck et alj l200ll) ^. The theoretical model 
is applicable out to halo masses of Mvir < W^'^Mq, giving rise 



to a cut-off in the predicted c. 



■2 X lo-i 



-M. 



2xl0'* 



relation at M- 



2x10*' 



10 Mq. The best-fit power-law index is a ~ 0.2. The data appear 
to agree well with the theoretical model of .BuUock et al.. ( 2001.) . 



8 SUMMARY 

The shape of the mass distribution in virialized halos provides a 
sensitive test of the CDM paradigm. In this paper we have analysed 
the mass profiles of a sample of galaxy clusters. The main results 
are as follows: 

• The mass distribution in the majority of the objec ts studied 
are consistent with the NFW model ( Navarro et al. '199^. Four ob- 
jects in the sample (2A 0335+096, Abell 478, PKS 0745-191 and 
Sersic 159—03) exhibit a flatter core. 

• The gas mass fraction measured in the central region of the 
clusters with flat core density profiles is high, indicating that the 
total mass may be underestimated within the central ~ 30 kpc 
of these objects. If this is the case, then either a significant non- 
thermal pressure component exists in the core, or the X-ray emis- 
sion has been incorrectly modelled. For 2A 0335+096 the central 

® Substituting r' = r/r^ into Equation 171 gives M(r')/M^-^\^ = 
[ln(l + CAr') - CAr7(l + c^r')] / [ln(l + ca) - ca/(1 + ca)]. 
^ The model data are found using the Fortran codes provided on J. Bul- 
lock's webpage. The theoretical Cvir-A/vir relation is scaled to C2x]^g4— 
^2x10" "=*ing *e definition M2X104 = M^ir x /(c2x lo* )//(cvir). 
where /(c) = ln(l + c) — c/{l + c), Cgxio^ = Cvir x 

(AvirA/2xl04/^2xl0lA^A^i,)^''^ and Avir = 104. 
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"•"-'-0.12 


6.6(4) 


Abell 3581 




(17+0-23 


9 81+6-30 


"-°^-0.20 


1 ,, + 1.05 
'-'-'-0.79 


0.081H? 


"•"-0.3 


23+0-16 
"-^•'-0.06 


1.3(4) 


Abell 2029 




n <:n+0.94 


4 38+l-^'' 
+-^°-1.76 


9 C7 + I.45 
^-^^-0.59 


n ,,+0.22 
"--'-'-0.19 


19+0-02 
"1^-0.02 


,+2.7 
^•-'-2.3 


77+0.26 
"-'^-0.13 


2.9(5) 


RXJ 1532.9 


+3021 


0.81+";?f 

— U.oo 


-7-7-1-2.28 
^- "-2.28 


2 24+4.72 
^■^^-0.85 


"■'-'-0.14 


"-'■'-0.02 


2 9+1.3 
^•^-1.3 


"-^°-0.15 


1.3(3) 


MS 2137.3- 


-2353 


34+0-74 


, 98+2.41 


1 74+1-22 
'-'^-0.42 


45+0-35 
"-^•'-0.30 


(, ,5+0.02 
"l-'-O.Ol 


5 0+1-^ 
-'•"-1.4 


52+0-22 
"--"^-0.10 


0.4(2) 


Sersic 159- 


03 


"-^"-0.10 


c i<;+3.42 


, 77+0.70 
'-^^-0.30 


57+0-52 
"-■"-0.36 


11+0-01 

"-"-0.01 


l-''-0.5 


,c+0.14 
"-^^^-0.07 


2.5(3) 



Table 4. Best-fit NFW model parameters for A = 200 and A = 2 X 10^. Also shown are the scaled radii enclosing a density contrast of A = 2500. 



- With Q a free parameter in the fits, a power-law density pro- 
file is ruled out at the 1 sigma level in Abell 1835 and Abell 
2029. (A power-law is also ruled out at the 1 sigma level in 2A 
0335+096, Abell 478, PKS 0745-191 and Sersic 159-03, al- 
though the inner slope in these objects is inconsistent with CDM 
predictions). Otherwise, we find no strong evidence for the ex- 
pected turn-over in the mass profile. For three objects in the sam- 
ple (Abell 3112, Abell 1795 and MS 2137.3-2353) the best-fit 
model is a power-law. 

- With a fixed at 1.0 the best-fit model is inconsistent with a 
power-law model for all objects except Abell 478. 

• The mass-concentration relation at a density contrast of 2 x 
10'' is cons istent with the theore tical prediction from A'^-body sim- 
ulations bv lBullock et alj j200ll) . 



9 DISCUSSION 

In the CDM scenario, simulations predict a universal form for 
the distribution of dark matter in virialized halos, independent of 
both m ass and cosmological parameters ilNavarro et alJl995Lll996L 
^93)- The universal density profile has a cuspy core, p oc r~", 
rolling over to p oc r ~^ at a characteristic radius known as the scale 
radius. The majority of the mass profiles obtained for the cluster 
sample in this paper are consistent with the NFW model; only four 
objects exhibit shallower density cores (2A 0335+096, Abell 478, 
PKS 0745 — 191 and Sersic 159—03). The discrepancy between the 
mass distribution measured in these systems and the NFW profile 
need not indicate any problem with the CDM paradigm, but may 
instead result from an underestimation of the core mass. The high 
gas mass fraction measured in the central region of all four systems 
with flat inner density profiles certainly suggests that this may be 
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Figure 10. Plot showing the dependence of concentrat ion on halo mass. The 
theoretical mass-concentration relation predicted bv lBullock et alJ <200lh 
(solid curve) and one-sigma confidence limits (dotted curves). Also shown 
is the best-fit power-law model to the data (dashed line). The colour-coding 
is the same as used in Figlsl 



temperature distribution is asymmetric, suggesting that the temper- 
ature and density obtained may not accurately represent the ambi- 
ent gas properties in the core for this object. 

• CDM simulations predict that the mass distribution in virial- 
ized halos differs from a power-law, turning-over at a characteristic 
radius, — rvir/c. 
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the case. Spuriously low mass measurements could be caused by 
the presence of non-thermal pressure support in the inner regions 
of these objects and/or incorrect spectral modelling of the emission 
from clusters with significant core substructure. Another possibil- 
ity concerns the importance of non-spherical effects on the mea- 
sured mass distribution. Havashi et al. 1 2004) have recently shown 
that rotation curves obtained from galaxies with triaxial halos may, 
in some cases, infer an erroneously flat core to the mass profile if 
spherical symmetry is assumed. The effects of non-axial symmetry 
on X-ray measurements of cluster masses have been investigated 
bv lPiffaretti et al. 12003.) . The authors show that the core mass will 
be underestimated if the halo is compressed along the line-of-sight. 
It is interesting that both Abell 478 and PKS 0745-191 both ap- 
pear highly elliptical. Analytic solutions to the hydrostatic equation 
for gas residing in a triaxial dark matter potential have been com- 
puted for isotherm al and polytropic temperature distributions by 
iLee & Sut j J2003h . These formula may in future work be applied 
to clusters to compute the mass distribution. 

Determining whether or not the flat cores measured in some 
clusters are real is of great importance to cosmological and clus- 
ter studies; if they are real, then either our understanding of dark 
matter is incorrect, or non-gravitational processes, such as dynam- 
ical friction acting on cluster galaxies moving through the core 
lEl-Zant et al. 2004), modify the mass distribution. 

The work carried out here has put some constraints on the 
mass distribution in a sample of galaxy clusters. Tighter constraints 
will be possible by imaging the clusters out to larger radii or, for the 
case of high redshift objects, by taking longer exposures, so that the 
radius at which the profile 'turns-over' may be robustly detected. 

In future work it will be important to investigate effects which 
might reconcile the observed flat density profiles found in this study 
with CDM model predictions. One possible cause f or the discrep- 
ancy is non-thermal pressure support in the core. ISanders et alJ 
i2005i) find inverse Compton emission from relativistic electrons 
which can contribute significant pressure near the centre of the 
Perseus cluster. Careful comparison of X-ray and lensing masses 
will also be important. 
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